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An algorithmic procedure for designing hybrid FIs/lJB digital filters 
is proposed and evaluated in this paper. The design is implemented as a 
two-stage optimization in which a Hooke and Jeeves optimization pro- 
cedure is used to optimize the iir component of the filter and the McClellan 
et al. optimization procedure is used to optimize the fir component of the 
filter. To evaluate this method, a set of eight low-pass filters were designed 
in which a single complex conjugate pole pair was used as the IIR com- 
ponent. The resulting designs were compared and contrasted with standard 
IIR low-pass filters and the optimal fir, linear-phase, low-pass filter in 
terms of multiplications, storage, and group delay properties. 

I. INTRODUCTION 

A wide variety of digital filter-design methods have been proposed 
and studied in the past several years. 1-9 Generally, these design methods 
can be classified as analytical or algorithmic, depending on the form 
of solution of the approximation problem which is used. Additionally, 
the resulting designs are classified as finite impulse response (fir) or 
infinite impulse response (im), depending on the filter properties, fir 
filters have the property that they can be easily designed to have 
exactly linear phase. Furthermore, linear-phase fir filters can be 
designed to approximate an arbitrary magnitude response to within 
given tolerances by using a sufficiently high order, iir filters cannot 
achieve a linear-phase response exactly, but are capable of approximat- 
ing sharp cutoff filters with considerably lower order filters than are 
required for the fir designs which meet identical magnitude specifica- 
tions. 10 Thus, for many practical filter applications, there is a trade-off 
between the exact linear-phase response obtainable using an fir filter 
and the reduced filter order obtainable using an iir filter. 

In this paper, a hybrid design is proposed which bridges the gap, 
somewhat, between the fir and iir filters. The hybrid filter is a par- 
ticular class of iir filter where the degree of the numerator of its system 
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Fig. 1 — Hybrid fir/iir filter. 



function is significantly higher than the degree of the denominator. 
Thus, it is reasonable to consider the hybrid design as a cascade of 
an Mth-order ira filter with Nth-order fir filter, as shown in Fig. 1, 
where N » M. In the designs discussed here, M is set to 2* and N is 
in the range of 27 to 51. The idea behind this hybrid design is to in- 
corporate the good features of both ira and fir filters, yielding a 
resulting filter which can meet arbitrary design specifications with a 
significantly smaller filter order than is required by the fir filter alone, 
and with a smaller group delay variation (i.e., more phase linearity) 
than is normally obtained by the iir filter alone. 

Based on the above discussion, the hybrid filter of Fig. 1 can be put 
in the form 

H(z) = 2W*)#rni(*), (!) 

where, by assuming a second-order denominator as mentioned above, 



and 



H U r(z) - (1 _ pe je z -i K i _ pe -^-i) 



Hfir(z) = £ h F (n)z- n , 

n-0 



(2) 



(3) 



where p is the radius of the pole in the z plane and 6 is the pole angle. 
(Although we have used a second-order iir filter in eq. (2) and in the 
examples of this paper, the design approach is general and can be 
applied to other orders of iir sections as well.) The overall filter is 
decomposed into a cascade of sections for two reasons. The first is to 
emphasize the fact that the hybrid filter is most readily realized as a 
cascade of an iir filter with an fir filter. In this way, the direct form 
structure can be used to realize the fir section so as to fully utilize 
the symmetry in the impulse response due to the linear phase condition. 

* Cases where M = 1 were also studied but have not led to any useful results. 
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The second reason that this cascade of sections is used is because it 
separates the design problem into two distinct parts; one which is 
fairly simple and one which can be solved using a well-known fir 
approximation method. 11 

To obtain the best second-order denominator hybrid filter approxi- 
mation to the desired specifications, the parameters p and of the iir 
filter are systematically varied, and an optimal fir filter is obtained 
for each set of parameters using the design program of McClellan et al. 11 
An alternative method of obtaining such a hybrid design would be to 
use a different type of algorithmic procedure (e.g., Ref. 8 or 9) and 
simultaneously vary both the numerator coefficients and the denomina- 
tor coefficients. Unfortunately, since the numerator order is so much 
higher than the denominator order, these optimization methods are 
not always successful. 

In the remainder of this paper we discuss the hybrid filter-design 
algorithm and present some typical results on low-pass filter designs 
obtained with this method. Then we compare and contrast the resulting 
designs with some standard iir low-pass filters including Butterworth, 
Chebyshev, and elliptic filters, and with the optimal fir filter. The 
bases of comparison are the number of multiplications per sample, the 
group delay variation, and the storage requirements. We conclude 
with a discussion of the properties of the hybrid filters. 

II. DESIGN ALGORITHM 

Since the algorithm is being applied to the design of low-pass filters 
(although it is equally useful for any arbitrary magnitude function), 
the desired frequency response of the system Hi(e ju ) is of the form 

where F p and F a are the passband and stopband edges, respectively. 
The tolerances are 8 P in the passband and 8„ in the stopband. Thus, 
|-ffr(e**)|, the magnitude response of the composite system of Fig. 1, 
satisfies the inequalities 

1 - 8 P ^ \Hi(e*) | ^ 1 + 8 P ^ a> ^ 2tF p 

£ |ff/(e*)| ^5. 2tF. £«£*-. (; 

The method in which the individual magnitude responses of the iir 
filters and the fir filters are chosen is shown in Fig. 2. Initial values 
are chosen for the iir coefficients p and 9, based on the location of the 
highest Q-pole pair of an elliptic filter which meets the tolerance scheme 
of eq. (5). The reason for this choice will be discussed later. 
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Fig. 2 — Block diagram of the design algorithm. 

Based on the initial values of p and 0, a set of design specifications 
for the fir filter is obtained from eqs. (1) and (5) of the form 



1 - 5 



1 + 5, 
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(6) 



Equation (6) specifies the appropriate parameters (the band edges, 
desired values, and weighting) for the McClellan et al. arbitrary magni- 
tude fir design program. 11 The final parameter required is N, the 
filter length. An initial guess of the value of N is made, and the design 
program yields the optimal fir filter for the given specifications. A 
control loop is used to find the smallest value of N which can be used 
to meet the given input specifications. 
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Once the minimum 7V (for the given values of p and 6) is obtained, an 
outer loop optimization program is used to vary p and 6 to obtain the 
minimum value of N as a function of p and 0. The algorithm used to 
find the optimum values of p and 6 is the well-known Hooke and Jeeves 
optimization. 12 Generally, the Hooke and Jeeves algorithm is capable 
of optimizing a continuous function of (several) continuous variables. 
However, for this problem the function [_N(p, 0)3 is not continuous, but 
instead is discrete (integer) valued. To handle the problems created by 
this discrete-valued function, fairly careful control over the variation 
of p and had to be maintained. To illustrate this point, Fig. 3 shows 
a typical contour of the variation of N as a function of p. 

The minimum of this function occurs at the point labeled A. How- 
ever, because N(p, 0) is discrete, it is possible for the optimization 
algorithm to prematurely terminate on flat regions such as those 
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Fig. 3 — Typical contour of JV as a function of p. 
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labeled B and C. To minimize these difficulties, a careful choice is 
required of the initial step sizes, Ap and A0, used by the Hooke and 
Jeeves algorithm to vary p and 9. For example, if Ap is too large, the 
algorithm may not find the "valley" of the function. Similarly, if Ap 
is too small, the algorithm may "hang up" on a flat region as it reduces 
its step size. Successful choices of Ap and A0 appear to be approxi- 
mately one-half the distance across the "valley" of the function (e.g., 
one-half the distance between the two points labeled D in Fig. 3). For 
the examples in the next section, the choices Ap = 0.005, Ad = 0.02 
(radians) yielded good results. Figure 4 shows another example of the 
variation of N as a function of p and 0. Here we can observe that there 
are two reciprocal minima of N as a function of p. If the optimum p 
outside the unit circle is found, the corresponding optimum value of p 
inside the unit circle can be found by taking the reciprocal of p . 

A factor that strongly affects the efficiency of the algorithm is the 
initial choice of p and 0. If accurate initial estimates of p and are used, 
they can result in a significant speed-up of the design method. As 
stated earlier, good estimates were found to be the location of the 
highest Q-pole pair of an elliptic filter which meets the same tolerance 
scheme. Other good estimates appear to be the highest Q-pole pair 
locations of the Chebyshev 2 and Chebyshev 1 designs which meet the 
same tolerance scheme. For greater assurance that an optimum has 
been found, several of these starting points may be tried to see if the 
algorithm converges to the same value of p and 0. 

Another factor that strongly affects the efficiency of the algorithm 
is the strategy of the control loop for varying JV (see Fig. 2). Generally, 
it was found that a good initial guess of N after changing p and is 
the previous value of N. An efficient strategy for increasing or de- 
creasing A r to find its new minimum then appeared to be a "tree" 
search (or log search). That is, N is incremented or decremented by an 
amount AJV, depending on whether the previous choice of N yields a 
design which meets the specifications in eq. (6). On the next trial, AiV 
is reduced by one-half, and the process is repeated until AN = 1. The 
search can be terminated sooner if at any stage the fir design is 
sufficiently close to the tolerance requirements (e.g., within 1 percent) 
given in (6). Other variations on this strategy are also possible. 

III. EXPERIMENTAL RESULTS 

Using the hybrid filter-design algorithm of Fig. 2, several low-pass 
filters were designed, ranging from narrow-band to wide-band designs. 
Table I gives the filter specifications (band edges and ripple tolerances) 
for eight low-pass filters, along with the resulting im pole position, the 
length of the fir section, and the length of an optimal linear-phase 
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Table I — Low-pass filter design 



Filter 
Example 


Filter Specifications 


Optimum Hybrid 
fir/iir Design 


Equivalent 
Optimum 
fir Design 


Highest Q-Pole 
of Equivalent 
Elliptic Design 


F P 


F. 


b„ 


S. 


P 


9(rad) 


N 


NfIR 


Pi 


e.(rad) 


1 
2 
3 
4 


0.05 
0.1 
0.15 
0.2 


0.075 
0.125 
0.175 
0.225 


0.01 
0.01 
0.01 
0.01 


0.001 
0.001 
0.001 
0.001 


0.9636 
0.9761 
0.9787 
0.9859 


0.3363 
0.6817 
0.9979 
1.3133 


39 
50 
51 
51 


109 
107 
106 
105 


0.9776 
0.9807 
0.9768 
0.9746 


0.3238 
0.6367 
0.9530 
1.2684 


5 
6 
7 
8 


0.05 
0.1 
0.15 
0.2 


0.125 
0.175 
0.225 
0.275 


0.005 
0.005 
0.005 
0.005 


0.0005 
0.0005 
0.0005 
0.0005 


0.9634 
0.9390 
0.9582 
0.9516 


0.3601 
0.6805 
0.9935 
1.3127 


31 
27 
31 
29 


44 
43 
43 
42 


0.9534 
0.9453 
0.9482 
0.9416 


0.3401 
0.6605 
0.9736 
1.2927 



fir filter which also meets the design specifications. The first four 
examples are filters with a fairly narrow transition band (0.025), 
whereas the last four examples are filters with a wider transition band 
(0.075). As seen in Table I, the reduction in length of the fir filter in 
the hybrid, from the optimal fir linear-phase filter, is on the order of 
1.5:1 to 2:1 — i.e., there is a fairly significant reduction in fir filter 
length. 

Figures 5 to 7 show examples of the magnitude responses of filters 
nos. 1, 5, and 2, respectively. Figure 5a shows the log magnitude 
response of the iir filter, Fig. 5b shows the log magnitude response of 
the fir section, and Figs. 5c and 5d show the linear and log magnitude 
responses respectively, of the composite filter. It can be seen from 
these figures that the response of the iir filter makes the requirements 
for the fir section significantly easier to obtain. For example, the 
required tolerance near a = t is on the order of — 15 dB in order that 
the composite response be more than 60 dB down at this frequency. 

Figures 6 and 7 show examples of some undesirable characteristics 
of the frequency response which can occur in the composite filter. In 
these cases, there is a ripple of the magnitude response in the transition 
band of the filter.* Since there is no real constraint on the composite 
filter response in the transition band, this behavior is not unreasonable. 
The question becomes one of whether or not a filter with a ripple in the 
transition band is acceptable. In general, the answer to this question 
is that it depends on the intended application. In some cases, this 
behavior is acceptable, in others it is not. 

Table II gives the approximate height of the transition band ripple 
for each example of Table I. In example 5 (see Fig. 6d), the transition 
band ripple is down by 17 dB and may be perfectly acceptable. In 
example 2 (see Fig. 7d), the ripple peak amplitude is about 1.06— i.e., 



* Note that the pure fir filter does not have such a peak in the don't care region. 
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Table II — Transition band ripple for Table I examples 



Example 


Height of Transition Band Ripple 


Magnitude 


dB 


1 
2 
3 
4 


1.06 
1.16 
1.52 


0.5 

1.29 

3.64 


5 
6 

7 
8 


0.14 
0.022 
0.11 
0.063 


-17 
-33 
-19 
-24 



it exceeds the passband maximum response by about 6 percent. Gen- 
erally, such large ripples render the filter useless in many applications. 
One way to combat this nonmonotonicity of the magnitude response 
in the transition band would be to constrain the pole of the iir section 
to lie within the passband of the filter. Another possibility would be to 
constrain the response of the fir section beyond the passband edge 
to guarantee monotonic behavior of the magnitude response. Either 
of these alternatives would lead to higher values of the fir filter dura- 
tion, thereby somewhat negating the gains of using the iir filter. 

IV. COMPARISON OF HYBRID FILTER DESIGNS TO OTHER 
CONVENTIONAL DESIGNS 

The hybrid fir/iir design represents a trade-off between an fir 
design and an iir design. The usefulness of such an approach depends 
strongly on how it compares with other conventional designs. To il- 
lustrate where in this scope of design alternatives a hybrid approach 
might be competitive, we have compared the hybrid examples designed 
in Section III to examples of five other conventional low-pass filter 
designs. These designs included the optimum finite impulse response 
(fir) design, the Butterworth (but) design, the Chebyshev type 1 
and 2 (chebi, cheb2) designs, and the Cauer elliptic (ellipt) design. 

In Fig. 8, the various designs are compared on the basis of the num- 
ber of multiplications required for their implementation. In the fir 
designs and the fir parts of the hybrid designs, the symmetry of the 
impulse response was exploited. In the iir designs, it was assumed that 
a conventional implementation of cascaded second-order sections was 
used and that coefficients of value 1 and 2 are not implemented with 
multiplies. Figure 8a shows the results for examples 1 to 4 in Table I 
and Fig. 8b shows the results for examples 5 to 8. 
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Fig. 8— Comparison of filter designs on the basis of the number of required 
multiplications. 

The results of the comparison in Fig. 8 indicate that the hybrid 
design is often more efficient than the Butterworth design and is 
slightly less efficient than a Chebyshev 2 design for the same magnitude 
specifications. 

A second criterion used for comparison was the group delay. Figure 9 
shows plots of the group delay for each filter design for example 8. The 
fir filter had the largest fixed delay at zero frequency, but its group 
delay was exactly flat (i.e., zero dispersion) across all frequencies. The 
hybrid design had a relatively flat component across most of the pass- 
band with most of its dispersion near the edge of the passband and 
within the transition band. The next least dispersive design appeared 
to be the Chebyshev 2 design, followed by the Butterworth, elliptic, 
and Chebyshev 1 designs. 
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Another measure of dispersion used to compare the designs was 
the maximum minus the minimum group delay across the passband. 
These results are plotted in Fig. 10. Obviously, the fir design had 
exactly zero dispersion in these examples. The next best contenders 
appeared to be the hybrid and the Chebyshev 2 designs, followed by 
the elliptic, Butterworth, and Chebyshev 1 designs. In three examples, 
the order of the Butterworth filter was larger than what could be ac- 
commodated by the available design program, so these results were not 
included. It was essentially the large required order of the Butterworth 
filters that prevented them from having favorable group delay 
characteristics. 

A final comparison was made on the number of data storage loca- 
tions (i.e., state variables) necessary for the implementation of the 
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Fig. 9— Comparison of group delays of various filter designs for example 8. 
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Fig. 10— Comparison of filter designs on the basis of group delay deviation (max- 
min.) in the passband. 



filters. These results are given in Fig. 11. In this respect, the fir 
designs and the hybrid designs did not compare as favorably to the 
recursive designs. 

By cross-comparing the above characteristics, it is possible to obtain 
a good insight into the various trade-offs involved in each filter design. 
Obviously, no single approach stands out as being superior over all 
other designs in all respects. The choice of a given design must be 
weighted according to the needs of a specific application. 
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Fig. 11 — Comparison of filter designs on the basis of required number of storage 
locations for state variables. 



V. CONCLUSIONS 

An algorithmic procedure has been proposed for designing hybrid 
fir/iir digital niters. The procedure is based on the use of a well- 
known fir design algorithm for designing the fir part of the filter, 
and it is coupled with a well-known optimization algorithm which is 
used to design the iir part of the filter. A set of low-pass filters designed 
in this way were found to have characteristics in between those of 
optimal fir designs and conventional iir designs. 
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Several properties of the design algorithm have been discussed in de- 
tail as well as methods of choosing initial starting points and techniques 
for speeding up the algorithm. Because of the nature of the fir/iir 
design, it was found that ripples could occur in the transition band. In 
some examples, these ripples were found to be objectionable, i.e., their 
amplitude exceeded that of the passband gain although the constraints 
imposed by the algorithm for the passband and stopband regions were 
completely met. Thus, an issue which needs further investigation is 
that of incorporating additional constraints in the algorithm to control 
the amplitude of the transition band ripples which can occur in such a 
hybrid filter design. 
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